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ABSTRACT 

We investigate non-axisymmetric low frequency modes of a rotating and magnetized 
neutron star, assuming that the star is threaded by a dipole magnetic field whose 
strength at the stellar surface, Bq, is less than ~ 10 12 G, and whose magnetic axis 
is aligned with the rotation axis. For modal analysis, we use a neutron star model 
composed of a fluid ocean, a solid crust, and a fluid core, where we treat the core as 
being non-magnetic assuming that the magnetic pressure is much smaller than the 
gas pressure in the core. For non-axisymmetric modes, spheroidal modes and toroidal 
modes are coupled in the presence of a magnetic field even for a non-rotating star. 
Here, we are interested in low frequency modes of a rotating and magnetized neutron 
star whose oscillation frequencies are similar to those of toroidal crust modes of low 
spherical harmonic degree and low radial order. For a magnetic field of Bq ~ 10 7 G, we 
find Alfven waves in the ocean have similar frequencies to the toroidal crust modes, and 
we find no r-modes confined in the ocean for this strength of the field. We calculate the 
toroidal crustal modes, the interfacial modes peaking at the crust/core interface, and 
the core inertial modes and r-modes, and all these modes are found to be insensitive 
to the magnetic field of strength Bo <10 12 G. We find the displacement vector of the 
core I' = \m\ r-modes have large amplitudes around the rotation axis at the stellar 
surface even in the presence of a surface magnetic field Bq ~ 10 10 G, where I' and m 
are the spherical harmonic degree and the azimuthal wave number of the r-modes, 
respectively. We suggest that millisecond X-ray variations of accretion powered X- 
ray millisecond pulsars can be used as a probe into the core r-modes destabilized by 
gravitational wave radiation. If the /' = \m\ = 2 r-mode is excited, we will have the 
pulsation of the frequency ~ 451/3 with Q being the spin frequency of the star. 
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1 INTRODUCTION 

Oscillation of strongly magnetized neutron stars has attracted an intense attention in recent years, particularly motivated 
by the discovery of quasi-periodic oscillations (QPOs) of magnetar candidates (e.g., Woods & Thompson 2006), which are 
believed to be one of the observational manifestations of global oscillations of neutron stars that have a strong global magnetic 
field of order of B ~ 10 15 G at the stellar surface (e.g., Duncun 1998, Israel et al 2005; Strohmayer & Watts 2005, 2006; 
Watts & Strohmayer 2006). Thus, recent theoretical studies of the oscillations of magnetized neutron stars have been mainly 
concerned with the stars having extremely strong surface magnetic fields So >10 G, and these studies have suggested that 
the QPOs observed in the magnetar candidates are attributable to the toroidal crust modes and Alfven modes of the neutron 
stars (e.g., Piro 2005; Glampedakis, Samuelsson & Andersson 2006; Lee 2007, 2008; Sotani, Kokkotas & Stergioulas 2008; 
Cerda-Duran, Stergioulas & Font 2009; Colaiuda, Beyer & Kokkotas 2009; Bastrukov et al 2009; Sotani & Kokkotas 2009). 
For toroidal modes of strongly magnetized neutron stars, however, it is intriguing, from the theoretical point of view, that Lee 
(2008) obtained discrete frequency spectra of the magnetic modes, but Sotani, Kokkotas & Stergioulas (2008), Cerda-Duran, 
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Stergioulas & Font (2009), and Colaiuda, Beyer & Kokkotas (2009), using a different numerical method from that used by Lee 
(2008), suggested the existence of continuum frequency spectra of the modes, as originally discussed by Levin (2006, 2007). 

The burst oscillation observed in low mass X-ray binaries (LMXBs) can be another example in which a magnetic field 
plays an important role in the oscillations of neutron stars, although the strength of the field at the surface of the neutron 
stars in LMXBs is thought to be less than ~ 10 10 G, much weaker than that for the magnetar candidates. For the burst 
oscillation, the hot spot model (e.g., Strohmayer et al 1997; Cumming & Bildsten 2001; Cumming et al 2002) and the Rossby 
wave model (Heyl 2004; Lee 2004) have been proposed, but it seems none of the models is accepted as the one that fully 
explains the observational properties of the burst oscillation. In the Rossby wave model, Heyl (2004, 2005) , Lee (2004), and 
Lee & Strohmayer (2005) have examined the possibility that the burst oscillation is produced by low frequency Rossby waves, 
called r- modes in the astrophysical literature, propagating in the surface fluid region (ocean), disregarding the effects of a 
magnetic field on the low frequency waves. We note that Bildsten & Cutler (1995) discussed the modal properties of g modes 
propagating in the fluid ocean above the solid crust of mass accreting neutron stars as a possible mechanism responsible for 
the ~ 6Hz QPOs observed in LMXBs. In their paper, on the assumption that the magnetic pressure, pb, is much smaller 
than the gas pressure, p, the critical strength of a magnetic field, B c , below which the magnetic field has no significant effects 
on the (j- modes, was estimated to be B c ~ 2 x 10 10 G for the accreted envelopes composed of carbon (see also Piro & Bildsten 
2005). Since their argument is based on a local analysis and on the assumption ps <S p, we think it useful to carry out a 
global modal analysis of low frequency waves propagating in the magnetized fluid ocean of a neutron star. 

Accretion powered millisecond X-ray pulsars in LMXBs show small amplitude, almost sinusoidal X-ray time variations, 
the dominant period of which is thought to be the spin period of the underlying neutron stars. Lamb et al (2009) argued 
that the millisecond X-ray variations are produced by an X-ray emitting hot spot located at a magnetic pole rotating with 
the star, and that so long as the symmetry center of the hot spot is only slightly off the rotation axis the X-ray variations 
produced have small amplitudes and become almost sinusoidal. They also suggested that if the hot spot is located close to 
the rotation axis, a slight drift of the hot spot away from the rotation axis leads to appreciable changes in the amplitudes and 
phases of the X-ray variations. Lamb et al (2009) pointed out that a temporal change in mass accretion rates and hence the 
radius of the magnetosphere, for example, can cause such a drift of the hot spot. We think it is also interesting to consider 
the possibility that global oscillations of the neutron stars work as a mechanism that perturbs the hot spot periodically. 

In this paper, to calculate global oscillations of a rotating and magnetized neutron star, we use the method of series 
expansion, in which the angular dependence of the perturbations is represented by series expansion in terms of spherical 
harmonic functions with different spherical harmonic degrees I for a given azimuthal wave number m (e.g., Lee & Strohmayer 
1996, Lee 2007). We calculate low m, low frequency modes of a neutron star composed of a fluid ocean, a solid crust, and a 
fluid core, where the star is assumed to be threaded by a dipole magnetic field, the strength Bo of which at the surface is 
smaller than ~ 10 12 G. The method of calculation we employ is presented in §2, and the numerical results are described in 
§3, and we discuss a local analysis for low frequency modes in the magnetized fluid ocean in §4, and we give discussion and 
conclusion in §5. 



2 METHOD OF CALCULATION 

We consider small amplitude oscillations of rotating and magnetized neutron stars in Newtonian dynamics, and no general 
relativistic effects on the oscillations are considered. We employ spherical polar coordinates (r,9,<j>), whose origin is at the 
center of the star and the axis of rotation is given by 6 = 0. We assume a dipole magnetic field given by 

B — p m X7(cos9/r 2 ), (1) 

where p, m is the magnetic dipole moment. For simplicity, we also assume that the magnetic axis coincides with the rotation 
axis. Since the dipole field is a force-free field such that (V x B) x B — 0, the field does not influence the equilibrium structure 
of the star. When we assume the axis of rotation is also the magnetic axis, the temporal and azimuthal angular dependence of 
perturbation can be represented by a sing le factor e i(m '*+ a ' t) , where m is the azimuthal wavenumber around the rotation axis 
and lo = a + m!) is the oscillation frequency observed in the corotating frame of the star, where a is the oscillation frequency 
in an inertial frame and fi is the angular frequency of rotation. The linearized basic equations used in a solid region of the 
star are given by 

1 n 1 

- uj 2 (, + 2iujSl x £ = -V • <r' - %V • a + - — (V x B') x B, (2) 

P + V • (M) = 0, (3) 

<«> 

B' = V x (£ x B), (5) 
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where p is the mass density, p is the pressure, £ is the displacement vector, and p' , p' , and B' are the Euler perturbations of 
the density, the pressure and the magnetic field, respectively, and A is the Schwartzshild discriminant defined by 

_ dlnp 1 dlnp . . 

and 

<91np 



r 



dlnp 



(7) 

ad 

In equation (2), a' denotes the Eulerian perturbation of the stress tensor, which is derived by using the Lagrangian perturbation 
of the stress tensor defined, in Cartesian coordinates, by 

Scjij = (Tipu) Sij + 2p (uij - ^uSij^J , (8) 
where is the strain tensor given by 

l (db , dij 



2\dx j + dx i )' (9) 

and n is the shear modulus and u = Ef =1 itn, and Sij is the Kronecker delta. Note that we have employed the Cowling 
approximation neglecting the Eulerian perturbation of the gravitational potential, and that no effects of rotational deformation 
are included. 

We can obtain the equation of motion for a fluid region by replacing the terms <x and a' in equation (2) by — p5,j and 
—p'Sij, respectively, and we do not need to consider equations (8) and (9) for fluid regions. 

Since the angular dependence of perturbations in a rotating and magnetized star cannot be represented by a single 
spherical harmonic function, we expand the perturbed quantities in terms of spherical harmonic functions Y™ (9, (f>) with 
different Is for a given m, assuming that the axis of rotation coincides with that of the magnetic field. The displacement 
vector £ and the perturbed magnetic field B' are then represented by a finite series expansion of length j m ax as 

Jmax 

* = 2 [ S h (r)Y™(e, 4>)e r + H h (r)VuY^(e, <f>) + T,, (r) e r x VhY^, 0)J , (10) 



3 = 1 

and 
B' 



Imax 

[b^(r)Y l T(e,4>)e r + bP j (r)VnY l T(e, C j>) + bT j (r) e r x V H F™(^)] e 1 "', (11) 



B (r) 

3 = 1 

and the pressure perturbation, p' , for example, is given by 

p' = E^M^(^K" (12) 

3 = 1 

where 

Vh = ee iL + e <t>^—z-ii> ( 13 ) 
o6 sin o<j> 

and Bo(r) — p, m /r 3 , and lj = \m\ + 2(j — 1) and l'j — lj + 1 for even modes, and lj = \m\ + 2j — 1 and l'j = l 3 ; — 1 for odd modes, 
respectively, and j = 1, 2, 3, • • • , j m ax- Substituting the expansions (10), (11) and (12) into the linearized basic equations 
(2)~(5), and (8), we obtain a finite set of coupled linear ordinary differential equations for the expansion coefficients . , , 
T;/ and so on, which we call the oscillation equation. For non-axisymmetric modes with m / 0, we cannot expect decoupling 

3 

between spheroidal (polar) modes and toroidal (axial) modes even for fi = 0, although for axisymmetric modes with m = 0, 
spheroidal modes and toroidal modes are decoupled when Q — (e.g., Lee 2007). 

In this paper, we treat for simplicity the fluid core being non-magnetic, which may be justified if we assume the magnetic 
pressure is much smaller than the gas pressure in the core. The oscillation equations used in the magnetized regions are given 
in the Appendix, in which the jump conditions imposed at the solid/fluid interfaces and the boundary conditions applied at 
the stellar center and the surface are also given. The oscillation equation used in the non-magnetic fluid core is found, for 
example, in the Appendix of Lee & Saio (1990). 



3 NUMERICAL RESULTS 



For the neutron star model used for modal analysis, we employ a cooling evolution model called NS05T7, which was calculated 
by Richardson et al (1982). The mass of the model is M ~ O.5M and the central temperature is T ~ 10 7 K, and the model 
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Figure 1. Low frequency \m\ = 1 modes of the model NS05T7 arc plotted as functions of — sgn(m)f2 for the case of Bo = 10 7 G for even 
modes in the left panel and odd modes in the right panel, where sgn(m) = m/\m\ and the positive and negative sides of the horizontal 
axis are for prograde and retorgrade modes, respectively. 



is composed of a fluid ocean, a solid crust and a fluid core, and because of density stratification (/-modes propagate in the 
fluid regions. The more detailed properties of the model, such as the equations of state used, are described, for example, in 
McDermott, Van Horn & Hansen (1988), who carried out modal analysis of the model assuming no rotation and no magnetic 
fields. Since this is a low mass model, the solid crust is rather thick, and the low radial order crustal modes of low spherical 
harmonic degree I are well separated from the / and p modes of the model, which makes the modal analysis much simpler 
than the case in which the crustal modes have frequencies similar to those of the / and p modes as expected for more massive 
neutron stars with accreted envelopes. 

Let us begin with the case of a weak magnetic field of strength Bo = 10 7 G. In the presence of a magnetic field, the 
existence of Alfven waves, whose oscillation frequency depends on the strength of the magnetic field and on the direction of 
wave propagation relative to the magnetic field, inevitably makes frequency spectra complex. For the case of B = 10 7 G, the 
frequency spectra of low frequency modes are too complicated to calculate completely, as suggested by Figure 1, which plots 
the oscillation frequency u> of low frequency \m\ = 1 modes as functions of Q for even modes (left panel) and odd modes 
(right panel), respectively, where we have used j m ax = 10. We notice that there appear in Figure 1 several kinds of modes, 
that is, the inertial modes and r-modes in the fluid core, the toroidal crust modes, the interfacial modes whose £ r amplitudes 
peak at the interface between the core and the solid crust, and the modes confined in the fluid ocean. The inertial modes 
and r-modes are rationally induced modes, whose oscillation frequencies are approximately proportional to fi and are found 
on almost straight lines tending to the origin in the figure, where the r-modes appear on the side of retrograde modes. The 
loci of the oscillation frequencies ui as functions of fi for the core inertial modes, r-modes, the toroidal crust modes, and the 
interfacial modes are more clearly seen for the case of a stronger magnetic field of Bo = 10 G as shown by Figure 2 below. 
We note that the toroidal crust modes are insensitive to the magnetic field of strength Bo <10 12 (e.g., Lee 2007), and that 
the inertial modes and r-modes in the fluid core and the interfacial mode at the core/crust interface are not strongly affected 
by the magnetic field even if their eigenfunctions extend to the surface through the magnetized solid crust and fluid ocean. 

However, the modes that are confined in the surface ocean are strongly influenced by a magnetic field as weak as 
Bo ~ 10 7 G and show extremely complicated frequency spectra. Although there usually appear gravito-inertial modes and 
r-modes as non-axisymmetric low frequency modes confined in the ocean for a non-magnetized neutron star, we find no ocean 
r-modes in the presence of a magnetic field, even if it is as weak as Bo ~ 10 7 G. We also find that Alfven modes come into 
the frequency spectra of the low frequency modes, modifying the spectra of the gravito-inertial modes. We notice that there 
exist two branches of modes in Figure 1. In one branch of modes, the oscillation frequency u> increases as \fl\ increases, while 
the oscillation frequency decreases with increasing \Q\ in the other, although there occurs frequent avoided crossings between 
the two branches of modes as Q varies. Based on the local analysis presented in §4, we think the former can be regarded as 
(gravito-)inertial modes and the latter as Alfven modes. However, we have to note that the eigenfrequencies of the modes in 



Non-axisymmetric low frequency oscillations of rotating and magnetized neutron stars 5 



o.i 



0.08 



0.06 



0.04 



0.02 



t 

l'=4 n = 



I 

1 = 1 2 



t 

l'=2 n=0 



-0.04 



-0.02 



0.02 



0.1 



0.08 ~ 



0.06 



0.04 



0.02 



t 

l'=3 n=0 



l"=1 n=0 



0.04 



-0.04 



-0.02 



0.02 



0.04 



-sgn (m) Si/ (QM/R 3 ) 172 



-sgn (m) Si/ (GM/R 3 ) 1/2 



Figure 2. Low frequency \m\ = 1 modes of the model NS05T7 plotted versus — sgn(m)f2 for the case of Bo = 10 10 G for even modes in 
the left panel and odd modes in the right panel, where sgn(m) = m/\m\ and the negative and positive sides of the horizontal axis are for 
retrograde and progradc modes, respectively. The symbols ;/t n , i«2, and ;/r n denote the toroidal crust modes, the interfacial mode whose 
amplitudes peak at the core/crust interface, and the core r-mode, respectively, where I and V denote the spherical harmonic degrees and 
the subscript n is the radial order, which is usually equal to the number of radial nodes of the dominant component of the eigenfunctions. 



the two branches do not necessarily reach good convergence even if j max is increased to j max ~ 20. In this sense we are not 
sure that the modes confined in the surface ocean we calculate are discrete modes with real frequencies. 

In Figure 2, we plot low frequency \m\ = 1 modes against £1 for the case of B = 10 10 Gfor even modes (left panel) and odd 
modes (right panel), where we have used j max = 10. In the figure, the symbols ;/£„, ii'2, and ;/r n respectively denote the toroidal 
crust modes, the interfacial mode whose £ r amplitudes peak at the core/crust interface, and the core r- modes whose frequency 
lu tends to 2m£l/ [I' (I' + 1)] as — > 0, where I and I' denote the spherical harmonic degrees and the subscript n indicates 
the radial order, usually corresponding to the number of nodes of the dominant eigenfunction in the propagation region. Note 
that we do not attach any labels to inertial modes in the fluid core whose frequencies are approximately proportional to O 
and are found on almost straight lines tending to the origin of the figure. Because the fluid core is almost isentropic such that 
the Brunt- Vaisala, frequency N is extremely small, the odd r-mode labeled ii = \ m \r n =o is the only r-mode we can find for a 
given value of m (e.g., Yoshida & Lee 2000a). We note that the eigenfrequencies of the modes plotted in Figure 2 obtain good 
convergence when j max is increased to j max >10. We have calculated low frequency \m\ = 1 modes for B — 10 12 G in the 
same frequency range as that in Figure 2, and obtained the result quite similar to that for the case of Bo = 10 10 G. 

Because of the effects of rotation, the frequencies w of the toroidal crust modes and the interfacial mode vary as £1 
changes. If we expand the oscillation frequency of a mode as lu = u>o + mCiQ. + 0(Q 2 ), where loq is the frequency of the mode 
for fl = 0, the Ci coefficient describes the first order response of the oscillation frequency to small Q, and a list of the C\ 
coefficients of rn — 2 spheroidal modes of the model NS05T7 is tabulated in Lee & Strohmayer (1996), where no magnetic 
effects are considered. For the toroidal crustal modes, on the other hand, Lee & Strohmayer (1996) showed Ci = 1/ [V (I' + 1)] 
for non-magnetized stars. Note that for rotationally induced modes such as inertial modes and r-modes, we have loo = 0. 
Figure 2 indicates that the frequency behavior of the crust modes for slow rotation is consistent to that expected from the 
Ci coefficients, even in the presence of a magnetic field of strength Bo ~ 10 10 G, although it is also clear that as increases 
from f2 = 0, the deviation of to from the expansion loo + mCifl quickly becomes significant. The deviation of the frequency 
from the expansion may be partly caused by avoided crossings between two different modes. An example found in Figure 2 is 
the avoided crossing between ; / =2 i„=o and 1=1X2, which is a crossing between a toroidal mode and a spheroidal mode. 

In Figure 3, we plot the eigenfunctions xSi 1 (x), xHi 1 (x), and xiT^(x) of several low frequency \m\ = 1 modes for the 

case of Bo = 10 10 G and tt/^GM/R 3 = 0.04, where x = r/R, and the low frequency modes we plot are the even toroidal 
crust mode j' = 2in=o, the even interfacial mode i=\ii, the odd r-mode ;'=ir n =o, and the odd toroidal crustal mode ;/ =3 t rl= o , 
and except for the r-mode, which is a retrograde mode with m = 1, all the other modes are prograde modes with m = — 1. 
Note that the toroidal crustal modes have appreciable amplitudes in the fluid core due to the effects of rotation. We find the 
eigenfunctions of the modes in the presence of the dipole magnetic field of Bo = 10 10 G are quite similar to those found in the 
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Figure 3. Eigenfunctions xSi 1 (panel a), xHi 1 (panel b), and xiT^ (panel c) plotted against x = r/R for the low frequency |m| = 1 
modes for the case of B = 10 10 G and Sl/^GM/B? = 0.04, where WC USG J max 

= 10. In each panel, the solid line is for the even toroidal 
crustal mode ; / =2 *n=0, the dashed line for the even interfacial mode ; = ii2, the dash-dotted line for the odd r-modc ;/ = 1 r n =o> and the 
dotted line for the odd toroidal crustal mode ('=3*^=0- Note that the amplitude normalization is given by xS^ = 1 at the stellar surface. 
Except for the r-mode, which is a retrograde mode with m = 1, all the other modes are prograde modes with m = — 1. 



absence of the magnetic field (see, e.g., Lee & Strohmayer 1996, Yoshida & Lee 2001), suggesting that these modes are not 
strongly affected by magnetic fields of that strength. Note that in the vicinity of the stellar center, since the eigenfunctions 
Si t and Hi ± are approximately proportional to x ll ~ 2 (e.g., Unno et al 1989), the functions xSi 1 and xHi 1 behave as oc x° for 
even modes of h = \rn\ — 1 and as oc x 1 for odd modes of h — \m\ + 1 = 2, as indicated by the panels (a) and (b). 

From the observational point of view, it is useful to know the 9 dependence of the displacement vector of the modes at 
the surface of the star. As noted in the previous section, for rotating and magnetized stars, the eigenfunction of an oscillation 
mode cannot be represented by a single spherical harmonic function Y; m (9, <j>) and hence their surface pattern can be largely 
different from that for non-magnetic and non-rotating stars. To see the angular dependence of the displacement vector at the 
surface, we introduce the functions X J (8) defined by 

2max 

X r {0)j m += , £ i S lj {R)Y??{0,4>) (14) 
2=1 

Jmax 

X (9) e im * = e e -J2 [ H h (R)^nY^(9, 0) + T,, (R) e r x V H >7 (0, 0)] , (15) 
2=1 

X 4, {e) eim ^ = _ ie ^ . 'g {r)VhY ^ { ^ 0) + T ^ {R) ^ x VHy ™ (e> , 



2=1 



(16) 

In Figure 4, we plot the functions X j (6) of the low frequency modes presented in Figure 3 for the case of Bo = 10 10 G and 
\J GM/R 3 — 0.04, where in each of the panels the solid, dashed, and dotted lines in each panel respectively denote the 
functions X r , X 6 , and , which are normalized by their maximum amplitudes. Since the modes have low frequencies, the 
maximum amplitudes of £ r are much smaller than those of £e and that is, the horizontal and/or toroidal components of 
the displacement vector are dominant over the radial component. It is to be noted that although the amplitudes of the ocean 
r-modes tend to be confined to the equatorial regions (see, e.g., Lee 2004), this is not the case for the l' = \m\ r-modes in the 
fluid core, which have large amplitudes in the regions around the the poles. 

In Figure 5, we plot the frequencies lo of low frequency \m\ = 2 modes against f2 for the case of Bo = 10 10 G for even 
modes (left panel) and odd modes (right panel), where we have used j max = 10. For the case of m| = 2, the toroidal crustal 
mode ;/ =3 t n= o appears as an even mode, while ;/ =2 t n =o and ; / =4 t n =o as an odd mode. We note that no interfacial modes 
appear in the frequency range shown in the figure. Figure 6 shows the functions X 3 (9) versus cos 9 for the core I' = m = 2 
r-mode at SI/ y 7 GM/R 3 = 0.04 for Bo = 10 10 G (left panel) and Bo = 0G (right panel), where the functions, evaluated at the 
stellar surface, are normalized by their maximum amplitudes. This figure shows that the surface pattern generated by the 
core I' = m — 2 r-modc is not affected by the field as strong as Bo = 10 10 G, which is also the case for the core /' = m = 1 
r-mode. 
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Figure 4. Functions X r (6), X°{9), and X*(0) versus cos 9 for the low frequency \m\ 
0.04 for j max 



1 modes for the case of B = 10 10 G and 
= 0.04 for jmax — 10, where the functions are evaluated at the stellar surface, and the solid, dashed, and dotted lines 
denote the functions X r (0), X 9 (6), and X^(0), respectively. Panels (a), (b), (c), and (d) are for the even toroidal crust mode ;/_ 2 tn=o, 
the even interfacial mode (=ii2, the odd core ; / =1 r n =o mode, and the odd toroidal crust mode ;/ =3 t n= o, respectively, where only the 
r-mode is a retrograde mode with m = 1, and the others are prograde modes with m = — 1. The functions are normalized by their 
maximum amplitudes. The functions X r {9) and X*(0) of even (odd) modes arc symmetric (antisymmetric) with respect to the axis of 
cos0 = 0, while the function X s (9) of even (odd) modes is antisymmetric (symmetric). 



4 LOCAL ANALYSIS 

To gain an understanding of the low frequency modes confined in the shallow fluid ocean calculated for the case of Bo = 10 7 G, 
we apply a local analysis to low frequency modes in a fluid region of a rotating and magnetized neutron star. Note that a 
local analysis of waves in rotating stars with no magnetic fields is found in, for example, Unno et al (1989). In this section, we 
assume that the x and t dependence of the perturbed quantities is given by the function exp [i (fe • x + wt)], where k is the 
wavenumber vector. Using the equation of motion for a fluid region and equations (3) to (5), we can derive a set of equations 
used for the local analysis: 

- u 2 i + 2icun x £ = --ikp + ^rVp + -!- (ifc x B') x B, (17) 
p p l Aivp y ' 

p' + ip(k-$)+$- Vp = 0, (18) 

P ' = r lP ($-A + p'/p) , (19) 

B' = ifc x (£ x B) , (20) 



which are combined to give 
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Figure 5. Same as Figure 2 but for |m| = 2. 
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Figure 6. Functions X r (0), X e (0), and X*{6) at the surface versus cos 8 for the m, = 2 i> = 2r n= o mode for the cases of Bo = 10 10 G in 
the left panel and Bo = 0G in the right panel, where Cl/yJ GM/R? = 0.04 and j ma x = 10 are assumed. The solid, dashed, and dotted 
lines respectively denote the functions X r (9), X e {6), and X^(0), which are normalized by their maximum amplitudes. 



(k-B) 

4np 



(a +a A )k- — -B 

K ' 4-zrp 



(fc . €) + (*^CLB) fc . 



2kjf2 x £ - ia 2 (£ ■ A) fc - 



0, (21) 



where a = \[VvpJp is the adiabatic sound speed and a a = B/^f^np is the Alfven velocity and B = \B\. 

To make the local analysis tractable, we employ a local Cartesian coordinate system whose z-axis is along the radial 
direction, and we assume that ft — fl z e z neglecting the local horizontal component of the rotation vector. Equation (21) can 
be rewritten into a form W£ = with W being a matrix, and the condition det W = leads to the dispersion relation: 



w 6 + A 4 uj + A 2 uj + Aiui + A = 0, 



(22) 



where 
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Ai = - (a 2 + a\) k 2 - a 2 A k 2 cos 2 9 - (2Q Z ) 2 - f3 2 , (23) 



A 2 = (2a 2 + a 2 A ) a A k 4 cos 2 9 + a 2 k 2 H N 2 + (2Q, Z ) 2 



A\ = — 4a 2 A gk 2 k z Q z sin ip cos 6, 



■ i ■ i i k z B z k z 

a fc. + a A k cos # + — — 2 cos 8—- — 
k 2 B k 



a 2 A k 2 H + 2a 2 , fc 2 cosf?-^-^- + (2f2 z ) 2 
B k 



Aq = — a 2 a 2 A k 4 cos 2 9 I a^fc 2 cos 2 # + -prN z 1 — a^fc* cos^ 9-^-f3 z , 



k 2 



4,4 2 n a z 2 



B 2> 



(24) 
(25) 
(26) 



where g = GM r /r 2 , cos6> = (fc ■ B) / (kB), k = \k\, B H = y/ B% + B 2 , k H = y/k 2 + fc 2 , sinV = (fc x B) z / (k H B H ), and 
N — yJ—gA is the Brunt- Vaisala frequency, and 



P 2 = - d -^g = N 2 + 9 -. 

dz a 2 

Note that the term proportional to Ai breaks the symmetry given by lj (— ft z ) = ui (£l z ). 
If we assume B = and fl 0, the dispersion relation reduces to 

co 2 {co 4 - [a 2 k 2 + (2Q Z ) 2 + P 2 ] co 2 + [a 2 k 2 (2Q Z ) 2 + a 2 k 2 H N 2 + (2Q Z ) 2 /3 2 ] } = 0, 

the non-trivial solution of which is 

On the other hand, if we assume fl = and B =fc 0, the dispersion relation reduces to 
(a; 2 — a\k 2 cos 2 0) 



(27) 



(28) 



(29) 



l 2 ,2 2,2 , 2 2 , 2 | 2 , 2 2 a , Ar 2 \ . 2 , 2 f B z fcfj B% . 2 , \ ,2 

w - (a fc + a A k + (3 ) to + a k y a A k cos 6 + -^N j + a A k [ + "p-^J sm ^ ) P 



B 2 fc 2 B 2 



+^ fc4 S- ft </>cos 2 0/3 2 = 0. (30) 
fc 2 _B 2 



If we can further assume sinip = 0, the solutions of the dispersion relation are separated into 

2 2 i 2 2/i 

u> — a A k cos 9, 

corresponding to the Alfven waves, and to 



Lo 2 = ± <j (a 2 + ai) fc 2 + /3 2 ± J [(a 2 + a 2 ,) fc 2 + /3 2 ] 2 - 4a 2 fc 2 ( a\k 2 cos 2 9 + ^N 2 + 



fc 2 



B 2 



(31) 



(32) 



For the case of B / and fi / 0, it is difficult to analytically solve the dispersion relation (22) in general . Here, we 
numerically solve equation (22) , which can be rewritten, by normalizing various quantities, into 



_ 6 A±_ A A 2 _ 2 Ai_ A 



Q 2 



where 
A 4 



/3 2 



( ,a = - (P + 9) (^ fc ) 2 - 9 (^ fc ) 2 cos2 ^ - 4fi 2 - 



(33) 



(34) 



^ = (2p + g) q (Rk) 4 cos 2 61 + p (Rk) 2 ^N 2 + 4n 2 



-4i a 9 ru , v.i f=, k z k H B H . , , 
_ = _4g-(i?fc) n, TX _ sm ^co8< 



p (Rk) 2 ^ + q (Rk) 2 \ cn^O + T r: - 2cnstf 



fc 2 



S fc 



+ 



g (i?fc) 2 ^ + 2 9 (fifc) 2 cos 0-± y + 4H 2 



n 2 ' 



A 



7^r = -P9 (fife) cos 6> 



"0 

where 



q(Rk) 2 cos 2 9 + iV 2 



2/ nn 4 2 a B z [3 2 



q (Rk) cos 6* 



B 2 fi 2 ' 



(35) 
(36) 
(37) 
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a 2 a\ a - r2 N 2 GM = Q z _ co GM 

P= m^' ?= wf' W gs = ~v' z = ^' " = <v " = V-^- (38) 

To solve the dispersion relation for a given neutron star model with M and 7? and for a given magnetic field B, we need to 
supply with appropriate values the following parameters, 

- 2 k z B z 

p, q, N, n z , (Rk), cosd, siaip, -p (39) 

although we note a relation given by 

k z B z + k H -B H k z B z k H B H , ,. n , 

cos * = kB = TfiT + Tb- cos ^ (40) 

For a given value of cos we have 
cozi ._ ™0-^/k){B z /B) 

COSf - (k H /k)(B H /B) ' (U) 
and the parameters k z /k and B z /B must satisfy an inequality cos 2 tp < I, that is, 

(t) 2 + (t) 2 " 2cOS 4t^- C ° s2 ^ ^ 
which can be rewritten as 

4 + 4 < 1. (43) 
where 



a; = 



1 (k z B z \ 1 (k z B z \ 2 2 8 2 .20 

In the following discussions, instead of k z /k and B z /B, it will be convenient to use the parameters / and Of defined by 

x = xofcos9f and y = yof sin Of with < / < 1. (45) 

Using these parameters we have k z /k — / cos (Of — 0/2) and B z /B = / cos (Of + 0/2), and p, q, N 2 , Cl z , (Rk), cos 0, f, and 
0/ are the parameters we need to specify. 

For the fluid ocean of the model NS05T7, typical values of the parameters p and N 2 are found to be 

p~l(T s , iV 2 ~ 10 5 , (46) 

and a typical value of the parameter q, depending on Bo, is q ~ 1(T 6 for B ~ 10 7 G and g ~ 1 for B ~ 10 10 G. Since 
ku/k — y/l — f 2 cos (Of — 0/2) <JC 1 for low frequency modes, we need / ~ 1 and 6/ ~ 6</2. In the following discussions, for 
simplicity, we assume / = 0.99 and Of = 0/2, and g — gs- Examples of numerical solutions of the dispersion relation (22) are 
given in Figure 7, where uo versus Cl z is plotted for cos# = 0.1 (dotted line), 0.5 (solid line), and 0.9 (dashed line) in the left 
panel, and lu versus €l z for q = 10~ 8 (dotted line), 10~ 6 (solid line), and 10 -4 (dashed line) in the right panel. It is interesting 
to note that the asymmetry due to the term proportional to Ai is too weak to become noticeable in the figure. As shown by 
the left panel, for a given value of cos#, the solution have two branches in this frequency region, and the upper and lower 
branches respectively correspond to gravito-inertial waves, for which to oc 2\Q, Z \ when \ fl z \ becomes large, and Alfven waves, 
for which uo oc a Ak cos 0. The minimum frequency in the inertial mode branch and the maximum frequency in the Alfven 
mode branch, which occur at Q z ~ 0, increase as cos# increases. These two branches of modes are reminiscent of the low 
frequency waves in the ocean plotted in Figure 1. We think that the two different mode branches associated with a mean 
angle between B and k appear as a pair of mode branches in Figure 1 and that the minimum and maximum frequencies at 
Q/tJGM/R 3 ~ in the pair depend on this mean angle, which may vary from one pair to another. As shown by the right 
panel of Figure 7, for a given value of cos 0, the frequency of the Alfven modes increases as q increases, and the inertial branch 
tends to the relation given by u> = 2\Q Z \ in the limit of q — > 0. 



5 DISCUSSION AND CONCLUSION 

Lamb et al (2009) proposed that the small amplitude, almost sinusoidal millisecond X-ray pulsation observed in accretion 
powered millisecond X-ray pulsars may be well explained by the hot spot model, in which the hot spots are assumed to be 
located at the magnetic poles, which are nearly aligned with the rotation axis. As discussed by Lamb et al (2009), even a small 
drift of the hot spot could produce appreciable changes in pulsation amplitudes of the X-ray pulsation. If this proposition is 
correct, it is interesting to pursue a possibility of using the millisecond X-ray pulsation to probe the core r- modes excited 
by gravitational wave radiation (Andersson 1998; Friedman & Morsink 1998). In fact, if the core r-modes with \m\ > 2 are 
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Q / (GM/R 3 ) 1/2 Q. / (GM/R 3 ) 1/2 

z z 

Figure 7. Left panel: Low frequency solutions of the dispersion relation (22) plotted as functions of Q z for cos0 = 0.1 (dotted line), 0.5 
(solid line) and 0.9 (dashed line), where we have assumed p = 10~ 8 , q = W' 6 , N 2 = 10 5 , (Rk) = 10 2 , / = 0.99, f = 6/2, and g = g s .; 
Right panel: Low frequency solutions of the dispersion relation (22) plotted as functions of fi z for q = 10~ 4 (dashed line), 10~ 6 (solid 
line), and 10~ 8 (dotted line), where we have assumed p = 10~ 8 , cos6» = 0.5, N 2 = 10 5 , (Rk) = 10 2 , / = 0.99, 6 f = 6/2, and g = g s . 



excited by the emission of gravitational wave, since the I' = m = 2 r-mode, which is the most strongly destabilized mode 
among the r-modes (e.g., Lockitch & Friedman 1999; Yoshida & Lee 2000a), produces the surface displacement vector £ whose 
horizontal and toroidal components at the surface have large amplitudes around the rotation axis as shown by Figure 6, the 
hot spot could suffer periodic disturbance from the r-mode. Note that the r-mode induced temperature perturbation, the 
surface pattern of which may be proportional to X r (8) e 1 ™* 1 , might generate X-ray variations, the amplitudes of which should 
be very small. If we write the oscillation frequency of r-modes as 

lo/q = K + k 2 h 2 + o (o 4 ) , (47) 

the coefficient Ko for the modes is given by 

ko = 2m/ [l' (l' + l)] , (48) 

and the coefficient «2 may depend on the equation of state and the deviation from the isentropic stratification in the core 
(e.g., Yoshida & Lee 2000a,b), where O = O/Oo. Since the neutron star core is nearly isentropic such that iV 2 ~ 0, we only 
have to consider the I' = \rn\ r-modes, for which we have ui w Kofi = 20/ (\m\ + 1), and we obtain the frequency u> w 20/3 for 
m = 2 in the corotating frame of the star and the frequency a = u) — mQ. w —40/3 in an inertial frame. It may be interesting 
to point out that if the r-mode of I' = m = 1 is also excited by some mechanism, this r-mode can produce long period 
variations in an inertial frame since the inertial frame frequency a ~ for this mode. Although no detection of periodicities 
whose frequency is approximately equal to 40/3 has so far been reported, if we detect the periodicities produced by the core 
r-modes of /' = \m\ in the X-ray millisecond pulsation, we can use the frequency deviation given by AcD = id — kqQ. w K2O 3 
to derive information about the equations of state and the thermal stratification in the core. The detectability of the r-mode 
pulsation may depend on the thickness of the crust, the property of the fluid ocean, and the strength of the magnetic field 
and so on. 

We have calculated non-axisymmetric low frequency modes of a rotating and magnetized neutron star, where we used 
a neutron star model composed of a surface fluid ocean, a solid crust, and a fluid core. We have assumed that the star is 
threaded by a dipole magnetic field but the fluid core can be treated as a non-magnetic region. For this model, we found that 
for a magnetic field of strength Bo ~ 10 7 G, Alfven waves in the surface ocean come in as low frequency modes, which largely 
modify the gravito-inertial modes in the ocean. However, the oscillation frequencies of the modes in the ocean are dependent 
on j'max, and they do not reach any good convergence even if ji max is increased to ~ 20. At present it is not clear that these 
ocean modes will converge to discrete modes with real frequencies in the limit of j max — » 00. We also found that no r-modes, 
which are confined to the surface ocean, can be found in the presence of the weak magnetic field. If this is also the case for the 
accreted envelopes expected in mass accreting neutron stars in binary systems, we need to reconsider the Rossby wave model 
for the burst oscillation in LMXBs (Heyl 2004, Lee 2004). On the other hand, the toroidal crust modes and the interfacial 
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modes at the core/crust interface, which show good convergence for j max >10, are found to be insensitive to magnetic fields 
of strength Bo <10 12 G. We also find that the core r-modes and inertial modes are not affected by the magnetic field even if 
their eigenfunctions extend to the surface through the magnetized crustal and envelope regions. However, this will not be the 
case for neutron stars having magnetic fields as strong as Bo ~ 10 15 G (e.g., Lee 2008). 

In the present paper, we employed for modal analysis a low mass neutron star model having a thick solid crust and a 
cold and thin surface ocean, for which the toroidal crust modes of low radial order and low spherical harmonic degree are 
well separated from the /- and p-modes. If we use more massive neutron stars with a hot accreted fluid envelope and a thin 
solid crust, the frequencies of the crust modes of low radial order and those of the /- and p-modes may overlap, leading to 
more complicated frequency spectra. We think it necessary to conduct similar modal analyses for such neutron star models 
to clarify the properties of the ocean modes in the presence of a magnetic field and to examine the possibility that the small 
amplitude millisecond X-ray pulsation can be used as a probe into the core r-modes. 



APPENDIX A: OSCILLATION EQUATIONS, JUMP CONDITIONS, AND BOUNDARY CONDITIONS 

In this Appendix, we present the oscillation equations solved for non-axisymmetric (m / 0) modes of rotating and magnetized 
neutron stars. As noted in the text (§2), we assume a dipole magnetic field whose axis is aligned with the rotation axis. 
We use Newtonian dynamics, and employ the Cowling approximation, neglecting the Euler perturbation of the gravitational 
potential. Using the dependent variables defined by 
p> 

(Vih=S h , {v 2 ) j = -J-, (y 3 ) 3 =H h , (y 4 ),=^., {b H )j=b$, (b T ),=6^ (6 S ) 3 =bf r , (Al) 
J grp J j 3 J o 

and 

2/ 5 = (Mib H +mA - 1 ib T ) /a, y 6 = (mK^b H + M () ib T ) /a, (A2) 

the oscillation equation for non-axisymmetric modes with m =fc for fluid regions threaded by a dipole magnetic field is given 
by 

^ = (^-3)2/ 1 -^3/ 2 + Ao2/3, (A3) 

d c lj 2 

r—r- = (rA + cjlu 2 ) y 1 + (1 - rA - U) y 2 - rnuciu) 2 y :i - vaij 2 Coy 4 + -^-r—R, (A4) 
dr x ' " 2 

M ( „1? + ™^ dj t = -\{k-¥^ + \r* y2 + ( M " - 5 Cl ) »" + - (A5) 

.-i dy., dy, IfV \ . _ x IV . _ x . _ x /„„ 1 _ \ 1 ib T 

" r ~dV + lT ~dT = 2 Vr7 " j Vl ~ 2 FT V2 + Vz + V " 2 C( V 2/4 " 2 a ~' (A6) 

dy. ( ._i MiC^X 2/2 (. 2MiAiM \ / 2m\ „ / dlna\ m ib T 

rJ iv = {-™ A ° + ^)^-^ + ( Lo — ) ^ - {" + -) M ^ + ( 2 -d^) + (A7) 

m\ ( 2m\ ,, /. 2m 2 . A / dlna\ 1 _ ib T , . „, 

a) K2/l - (" + —) M0V * + [ Ll - — Al J 2/4 + ( 2 - d^) V <> -2 Cl ~ (A8) 



dy (i ( m 
r—?- = ( v - 
dr 

with 



„ d (. b H ib T \ ( d\na\ ( b H ib T \ b s 

R = r— Co h rn - 2 Co h m - Co — (A9) 

dr y a a J V d In r / ^ a a J a 

and 

6 s = -AiK2/ 1 -2AiMo2/3-2m2/ 4 , (A10) 

where v = 2Q/u), and u) = w/ \J GM/R 3 with M and R being the mass and radius of the star, and 

dlnTW, dlnp (S/Rf a= c^pV Bg 

rflnr ' dlnr' 1 M r /M ' 4p B ' P 8tt ' 1 ' 

and Cq is the transposed matrix of Co- The non-zero elements of the matrices M , Mi, Co, Ci, K, A , and Ai that appear 
in the oscillation equation given above are defined by 



(Mok; = jjfi^+i. ( M o)i.i+i = i~^ J h+ 2 ' ( M 0« = ]Jji J h+i> (MiJi+i,,- = Jqr^7+2, (A12) 
(Co)j,j = + 2) J,™ +1 , (Co)j+i,j = (Ij + l)J™ + 2, (Ci)jj = ljJ™ + i, (Ci)j,j+i = + 3) J™+2, (A13) 
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(Kk< = 



L + V 



(K) 



J.3 + 1 — 



L + 2 



, (A ) w = Z,ft + 1) , (Ai)jj = + l) 



for even modes, and 



«i + 2 



(Co)i,i — ft - i)JiTi (Co)j,3+i — -ft + 2) jr+i, (Ci)jj - -ft + 1) jr, (Ci)j+i j — z^ jr+i, 



(K)w = — f 2 -, (KJi+i,,- = -fj-, (A k 3 = Z, ft + 1) , (Ai) w = Z; (Z; + 1) 



for odd modes, where 

J\. - 



I' 2 



4/ 2 - 1 

for I > \m\, and J ; m = for I < \m\. The matrices L and Li are given by 
Lo = 1 — mvA.Q 1 , Li = 1 — mi/A^ 1 



(A14) 

(A15) 
(A16) 

(A17) 

(A18) 
(A19) 



with 1 being the unit matrix. Note that lj — \m\ + 2(j — 1) and l'j = lj + 1 for even modes, and lj — \m\ + 2j — 1 and lj = l 3 ; — 1 
for odd modes. 

For a solid region threaded by the dipole magnetic field, on the other hand, we use the dependent variables defined as 



{zi) j =S lj , {z2)j=H h , (z 3 ) j =iT l ,, 
and 



+ 2a 1 ±(rz 1 ) + ^(Q ) b H +rmb T ), 

dr p v 7 



^ = (Ti - fai) J; (r 3 *i) - Ao* 2 

Z5 = Ql ( r ^l + Zl ) - ( Ml b H + rnA^ib T ) , z e> = ai r^ _ ^ (mA -i^ + 

\ dr J p y 7 dr p v 7 



where 

pi 2 4 

Ql = — , «2 = Ti — -Oi, «3 = Ti + -Qi. 
p 6 3 

The oscillation equation then becomes 

dzi 3ri a 2 1 2p B /_ , h . ., T\ 
r— — = zi H A ( )Z 2 H «4 CoO + mib , 

dr «3 Q3 «■> v 



Q3 Q?3P 



(A20) 

(A21) 
(A22) 

(A23) 

(A24) 



Mor ^l + mAr V^ = i fl + — ) Kzi + (Mo - ^Ci) z 2 + mA^zs ^-Kz 4 
dr dr 2 V «3 / V 2Q3 / 2Q3 



+ ( 

\a 3 p 



PB KCo - \l ) b H + -^-mK,V 



a 3 p 



(A25) 



i dz2 , dzz 
mA r— hMir— — 

dr dr 



i/, 3rA i 

-- H mA Zi + rn 

2 V CK3 / 



P B A- 1 /- J« 1 ft i 2PB 2 A -l\ .,T 

mAj Coo — - I 1 H m A ib , 

azp 2 \ a 3 p J 



(A26) 



dZ4 

dr 



(uV-4V-c 1 u, 2 V+^^)l 



Ct 3 



2p B 



CAiK 



Zl + 



mciuj 2 Vv\ + (v - 2ai - \ _ \ , M„ 



«3 / 



^2 



+ ( cioj Vv — m ) Coza 



\ Q3 / 



Z 4 + AqZ-s + 



4pb 



A() Mi - 2 f 1 - SH.) C 1 b * - ^ f 1 - ^ mi b T {A27) 
\ azJ J p V «3 / 



dr 



Z2 



(v- ^lEl) i + mcniV^Ao 1 + -^.MiAiK zi + -auj 2 VU + ^MiAiM,, - 2ail + 2ai ( 1 + — 1 A 

V as ) p p \ a a J 

8mPB ^ » - ^ 2 «A „ , o^ _ , 4 PB r„ , ^\,] b H + 4pB m / A -i _ £1 j\ . b T (A2g) 

V Q3 / 



p 



H- ! ■ - r.- + ^ ) Mi^3 - (l - ^L) z 4 + (V - 3) z 5 + 



4p B 
V 
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dz & 
dr 



-cioj Vv + 



4p B 



m I Kz\ + c\u> Vv + 



Sp B 



m M0Z2 + 



p V 2 / 

+ (V _ 3) * 6 + ^mAr 1 ^ + 1^ f M + i Cl ) l 6 I 
P p V 2 / 



(A29) 



We derive one set of the jump conditions imposed at the interface between the solid crust and the fluid ocean by assuming 
the continuity condition of the displacement vector at the interface (see e.g., Lee 2007): 



= 0, 



(A30) 



where [F = lim E ^o [F( r i + e ) — F{ r i — e)} . The other set of the jump conditions are derived from the continuity condition 

of the rr, r6, and r<j> components of the perturbed traction at the interface; 



[St,-j (r 



= 0, 



(A31) 



where Srij denotes the ij component of the perturbed traction, and the detail expression of 5rij is given in Lee (2007). As for 
the jump conditions at the interface between the solid crust and the fluid core, which is assumed non-magnetic, we use the 
continuity of the radial component of the displacement vector and of the rr, r6, and r<j> components of the perturbed traction 
at the interface. 

The surface boundary conditions for non-axisymmetric modes are given by 



Sp/p = 0, ib T = 0, b s + L + b" = 0, 



(A32) 

where 5p denotes the lagrangian perturbation of the pressure, and L + = 5ij + l) (Lee 2007). As for the inner boundary 
conditions at stellar center, we require that the functions ry 1 and ry 2 are regular at the center. 
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